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1.  INTRODUCTION 

The  disposal  of  obsolete  munitions,  propel¬ 
lants,  and  manufacturing  wastes  is  conducted  at 
Department  of  Defense  (DOD)  and  Department 
of  Energy  (DOE)  facilities.  The  most  common 
disposal  method  is  open  burning  (OB)  and  open 
detonation  (OD)  of  the  material,  which  occurs  in 
an  earthen  pit  or  bermed  area.  At  present,  the 
material  destroyed  in  a  single  detonation  typically 
ranges  from  100  to  5000  lbs,  whereas  the  quantity 
treated  in  a  bum  can  be  somewhat  larger  and  last 
from  minutes  to  an  hour.  OB / OD  activities  are  re¬ 
stricted  to  daytime  during  unstable  or  near-neutral 
atmospheric  stability. 

OB/OD  operations  generate  air  pollutants  and 
require  predictions  of  pollutant  concentrations. 
The  pollutants  include  SO2,  NOx,  particulates, 
volatile  organic  compounds  and  toxic  materials 
such  as  metals,  semivolatile  organics,  etc.  (An- 
drulis,  1992).  For  large  detonations  (1  -  3x  10“*  lbs), 
natural  dust  entrained  by  the  blast  is  an  additional 
contaminant.  Emissions  from  OB/OD  sources  have 
the  following  unique  features:  1)  “instantaneous” 
or  short-duration  releases  of  buoyant  material,  2)  a 
wide  variability  in  the  initial  cloud  size,  shape,  and 
height,  and  3)  ambient  exposure  times  from  clouds 
that  are  much  less  than  the  typical  averaging  times 
(~  1  hr)  of  air  quality  standards. 

Dispersion  models  are  used  to  estimate  pollu¬ 
tant  concentrations  given  the  source  and  meteo¬ 
rological  conditions.  However,  there  is  currently 
no  recommended  EPA  dispersion  model  to  address 
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OB/OD  sources.  The  most  widely-used  approach  is 
INPUFF  (Petersen,  1986),  a  Gaussian  puff  model, 
but  this  has  several  limitations  as  discussed  be¬ 
low.  Due  to  the  contraints  of  existing  models,  a 
model  development  program  was  initiated  under 
the  DOD/DOE  Strategic  Environmental  Research 
and  Development  Program. 

In  Section  2,  we  give  an  overview  of  the 
model  design  which  is  divided  into  “simple”  and 
“research”  components.  Sections  3  and  4  describe 
the  simple  component  which  includes  Gaussian  puff 
aind  analytic  plume  models.  This  development 
program  is  in  progress  and  is  currently  limited  to 
the  unstable  planetary  boundary  layer  (PBL). 

2.  MODEL  DESIGN  CONSIDERATIONS 


2.1  Background 

The  development  of  an  OB/OD  dispersion 
model  has  considered:  1)  the  limitations  of  existing 
models,  2)  current  knowledge  of  turbulence  aind  dis¬ 
persion  in  the  PBL,  and  3)  a  mobile  meteorological 
platform  under  development. 

Limitations  of  existing  models.  INPUFF  has 
been  used  to  model  OB/OD  sources  and  can  handle 
dispersion  from  individual  puffs  or  clouds  or  from 
a  sequence  of  puffs  in  a  short-duration  release.  Al¬ 
though  the  Gaussian  puff  approach  is  suitable  for 
OB/OD  sources,  INPUFF  has  the  following  limita¬ 
tions:  1)  It  adopts  dispersion  parameters  {ay,az) 
from  the  Pasquill-Gifford  (PG)  curves  or  from  Ir¬ 
win’s  (1983)  scheme.  2)  It  includes  Briggs’  (1971) 
plume  rise  expressions  which  apply  to  continuous 
releases  rather  than  to  instantaneous  sources  (puffs, 
clouds)  and  does  not  address  thermal  penetration 
of  elevated  inversions  capping  the  PBL.  3)  It  as¬ 
sumes  Gaussian  velocity  statistics  for  the  turbu- 


lence,  whereas  the  vertical  velocity  statistics  in  the 
unstable  PEL  are  positively  skewed  (Wyngaard, 
1988).  The  skewness  should  be  included  for  ver¬ 
tical  dispersion. 

For  OB/OD  sources,  the  PG  curves  are  defi¬ 
cient  in  that  they:  1)  are  based  on  dispersion  from 
a  ground-level  source  and  short  downwind  distances 
(<  1  km),  and  2)  are  selected  using  surface  meteo¬ 
rology,  which  does  not  account  for  the  PBL’s  ver¬ 
tical  structure.  For  large  detonations,  source  buoy¬ 
ancy  can  carry  emissions  to  several  100  m  or  the 
PBL  top;  one  must  then  deal  with  dispersion  over 
the  entire  PBL. 

PBL  turbulence.  Dispersion  in  the  PBL  de¬ 
pends  on  the  turbulence  length  and  velocity  scales 
which  differ  for  the  unstable  or  convective  boundary 
layer  (CBL)  and  the  stable  boundary  layer  (SBL). 
For  the  CBL,  the  length  and  velocity  scales  are  the 
CBL  depth  h  and  the  convective  velocity  scale  w,. 
Typical  values  of  m.  and  h  at  midday  over  land  are 
1-2  m/s  and  1-2  km.  Within  the  “mixed  layer" 
(O.l/i  <  z  <  h),  the  mean  wind  speed  and  turbu¬ 
lence  components — longitudinal  lateral  tuid 

vertical  a-uj — vary  little  with  height  z]  in  strong  con¬ 
vection,  (Tu,  (Tu,  ~  0.6w.. 

For  the  SBL,  the  turbulence  is  much  weaker 
with  eddy  sizes  proportional  to  z  near  the  surface 
and  typically  ~  10s  of  meters  or  less  in  the  upper 
part  of  the  SBL.  Models  and  observations  show 
that  the  velocity  scale  is  the  friction  velocity  u. 
(Wyngaard,  1988),  which  is  typically  ~  0.1  m/s  in 
strong  stable  stratification. 

Knowledge  of  the  PBL  turbulence  structure 
has  been  included  in  models  for  applications  (see 
Venkatram  and  Wyngaard,  1988). 

Mobile  meteorological  platform.  A  mobile  me¬ 
teorological  platform  is  being  developed  at  NOAA- 
ETL  to  obtain  the  PBL  variables  necessary  for 
modeling  since  many  DOD  facilities  are  in  remote 
locations.  The  platform  design  includes:  1)  a  radar 
wind  profiler  for  obtaining  the  three  wind  compo¬ 
nents  up  to  ~  3  km,  2)  a  radio  acoustic  sounding 
system  (RASS)  for  temperature  measurements,  3) 
a  mini-SODAR  for  measuring  winds  and  to  a 
height  of  ~  200  m,  4)  a  mini-lidar  system  for  ob¬ 
taining  the  PBL  depth  h,  and  5)  a  portable  meteo¬ 
rological  station  for  measuring  near-surface  winds, 
temperature,  turbulence,  and  heat  flux.  The  dis¬ 
persion  model  is  being  designed  for  efficient  use  of 
these  measurements. 

2.2  Overall  Model  Design 

A  model  h^archy  is  planned  including:  1)  a 
simple  computational  framework  for  routine  prob¬ 
lems,  and  2)  a  more  detmled  or  research  model  for 
nonroutine  problems.  In  the  simple  approach,  a 
Gaussian  puff  model  is  adopted  for  instantaneous 


sources  and  puff,  integrated-puff,  and  plume  models 
for  short-duration  releases.  For  the  research  frame¬ 
work,  a  Lagrangian  particle  and/or  puff  approach 
is  planned.  Both  frameworks  will  be  considered  for 
“onsite”  use  in  a  real-time  operational  mode  us¬ 
ing  data  from  the  mobile  meteorological  platform, 
i.e.,  for  day-to-day  decisions  on  OB/OD  operations. 
The  puff  and  plume  models  would  be  used  for  cli¬ 
matological  cuialyses  needed  in  risk  assessments. 

In  modeling,  the  important  aspects  to  address 
are:  1)  all  source-related  features  including  the  in¬ 
stantaneous  or  short-duration  nature  of  the  release, 
buoyancy-induced  rise  and  dispersion,  and  cloud  or 
plume  penetration  of  elevated  inversions,  2)  relative 
and  absolute  dispersion  expressions  that  explicitly 
include  PBL  turbulence  variables,  3)  meteorologi- 
cad  variables  including  their  vertical  profiles  from 
the  mobile  platform,  and  4)  a  treatment  for  puff 
and  plume  dispersion  about  complex  terrain. 

The  following  models  address  points  1  and  2 
above  and  must  be  expanded  to  include  points  3 
and  4.  Further  development  also  will  aiddress:  1)  a 
more  complete  description  of  initial  source  effects 
(detonation  cloud  size  and  height)  and  inversion 
penetration,  2)  a  more  complete  PBL  turbulence 
parameterization,  3)  averaging  time  effects  on 
concentration,  4)  the  entrained  dust  source  term, 
and  5)  deposition  of  gases  and  particles. 

3.  INSTANTANEOUS  SOURCES 


3.1  Dispersion  Model 


Concentration.  For  instantaneous  sources  or 
detonations,  a  Gaussian  puff  model  is  adopted  for 
the  short-term  mean  concentration  (C)  field: 


C  = 


Q 


{2T^Y/^arxCrry(^rz 


exp 


{x  -  Ut)"^  _  ip- 
2(7?,  2<t?, 


(z  -  hpp' 
2(7?.  . 


(1) 


where  Q  is  the  pollutant  mass  released,  U  is  the 
mean  wind  speed,  t  is  the  travel  time,  hg  is  the 
effective  puff  height,  and  cTrx,  <^ry,  and  Crz  are 
the  puff  standtu-d  deviations  or  relative  dispersion 
in  the  x,  y,  cind  z  directions,  respectively.  Here, 
hg  =  h,  +  Ah  where  h,  is  the  source  height  and 
Ah  is  the  cloud  rise  due  to  buoyancy;  x  and  y 
are  the  distances  in  the  mean  wind  and  crosswind 
directions. 

Currently,  we  are  considering  two  approaches 
for  estimating  the  peak  ground-level  concentration 
(GLC)  at  a  given  x:  1)  the  pealc  concentration  Cg 
in  the  elevated  buoyant  puff,  and  2)  a  pealc  found 
from  a  probability  distribution  of  concentration  at 
a  downwind  receptor.  The  Cg  is  the  puff  centroid 


rise  in  a  neutral  environment 


concentration  given  by  Cc  =  <5/[(27r)^/^CTriOrrj,nrj]) 
where  the  relative  dispersion  parameters  are  gener¬ 
ally  different  in  the  three  directions.  In  the  follow¬ 
ing,  we  assume  Oj-x  =  Ory  =  <^rz  —dr-  If  Cc  is 
used  as  an  estimate  of  the  peak  concentration,  an 
estimate  must  be  made  of  the  probability  of  it  be¬ 
ing  brought  to  the  surface;  one  possible  method  is 
given  by  Weil  et  al.  (1995). 

In  the  second  approach,  we  require  a  functional 
form  for  the  concentration  probability  distribution 
(e.g.,  a  gamma  distribution;  Deardorff  and  Willis, 
1988)  and  estimates  of  C  and  the  root-mean-square 
concentration  fluctuation  Oc  due  to  an  ensemble 
of  meandering  puffs.  The  probability  distribution 
and  the  Oc  model  remain  to  be  selected.  The 
C  field  including  puff  meandering  is  given  by  Eq. 
(1),  but  ndth  <7ra,  replaced  by  the  absolute 

dispersion  parameters — cTi,CTy,crz.  A  Gaussian 
distribution  for  C  is  applicable  to  the  SBL  where 
the  probability  density  function  (p.d.f.)  of  the 
vertical  velocity  w  is  Gaussian.  However,  for 
the  CBL,  a  skewed  w  p.d.f.  is  more  consistent 
with  laboratory  and  field  data.  A  skewed  p.d.f. 
is  adopted  here  and  is  parameterized  by  the 
superposition  of  two  Gaussian  distributions  (Weil, 
1988). 

The  C  field  due  to  an  ensemble  of  meander¬ 
ing  puffs  is  derived  from  pu,  following  the  same 
approach  as  applied  to  continuous  plumes  (Weil, 
1988).  The  resulting  expression  for  C  is 


_  (2) 
where  =  a^x/U  and  Zj  =  vjjxjU  with  j  =  1,2. 
The  Aj,  Oj,  and  wj  {j  =  1,2)  are  the  weight,  mean 
velocity,  and  standard  deviation  of  each  Gaussian 
p.d.f.  comprising  Pa,.  Equation  (2)  applies  for  short 
distances  such  that  the  plume  interaction  with  the 
ground  or  elevated  inversion  is  weak.  The  complete 
expression  for  C  includes  multiple  cloud  reflections 
at  the  ground  and  PBL  top. 

The  time-averaged  concentration  can  be  found 
from  the  dose  where  the  partial  dose  is  defined 
by  ■4){x,y,z,t)  =  f^C{x,y,z,t')dt'  and  the  total 
dose  by  V»oo  =  ''P{x,y,z,oo).  For  clouds  with 
short  passage  _times  over  a  receptor, _the  average 
concentration  C  can  be  obtained  from  C  =  — 

rp{ti))fTa,  where  the  averaging  time  Ta  =  h  —  ti.  If 
the  puff  passag^itime  AdrxIU  is  less  than  Ta,  then 

C  =  T/>oo/Ta. 

Cloud  rise  and  inversion  penetration.  Scorer 
(1978)  combined  theory  zmd  laboratory  experi¬ 
ments  to  obtain  the  following  expression  for  cloud 


Jlh  =  2.i5{MTt  + . 


(3) 


Mt  and  Ft  are  the  initial  momentum  and  buoy¬ 
ancy  of  the  cloud  and  are  given  by 


Mt  = 


and 


Ft  = 


qQt 

CpPa^a 


(4) 


where  Wo,  To,  and  Qt  are  the  initial  velocity, 
radius,  and  heat  content  of  the  thermal,  g  is  the 
gravitational  acceleration,  Cp  is  the  specific  heat  of 
air,  and  pa  and  Go  are  the  ambient  air  density  and 
potential  temperature. 

Scorer  also  found  the  puff  radius  to  be  r  = 
aAht,  where  A/i(  is  the  cloud  top  height  and  a 
is  an  empirical  entrainment  coefficient,  a  ranged 
from  0.14  to  0.5  with  a  mean  of  0.25.  The  relative 
dispersion  =  t/\/2. 

Using  field  observations,  Weil  (1982)  confirmed 
that  Eq.  (3)  was  a  good  fit  to  data  over  a  wide  range 
of  times.  Thus,  Eq.  (3)  is  suitable  for  the  initial 
rise  of  a  cloud,  i.e.,  before  it  is  limited  by  stable 
stratification.  The  Qt  can  be  determined  from 
the  mass  of  the  detonation  and  its  heat  content, 
H  =  1100  kcad/kg  TNT  equivalent. 

For  cloud  penetration  of  an  elevated  density 
jump,  results  have  been  found  from  laboratory  ex¬ 
periments  in  a  nonturbulent  environment.  Richards 
(1961)  obtained  an  empiricaJ  expression  for  the 
fraction  P  of  the  cloud  penetrating  the  jump;  P  ~ 
1  —  Q.bApi! ApTi,  where  Api  is  the  density  jump 
and  ApTi  is  the  average  density  excess  of  the  cloud 
when  it  reaches  the  jump.  The  ApTi  can  be  esti¬ 
mated  from  Ft  auid  r. 

The  Ap  in  a  detonation  cloud  is  related  to 
the  cloud  temperature  excess  A0  by  Apjpa  = 
AQ/Qa  with  A0  =  (3/47r)(5T/(PoCpr^)-  We 
can  then  rewrite  Richards’  expression  as  P  = 
1  —  {2ir/3){AQiPaCpa^h^ /Qt),  where  A0i  is  the 
temperature  jump  at  z  =  h.  The  h,  is  assumed  to 
be  zero  so  that  the  cloud  radius  at  the  inversion 
is  ah.  The  above  relationship  shows  the  strong 
sensitivity  of  P  to  ah. 

Figure  1  shows  examples  of  P  versus  the 
detonation  mass  W,  where  we  have  used  Qt  = 
W  •  H  and  a  =  0.25.  For  h  =  500  m,  one  can 
see  that  a  significant  fraction  of  the  cloud  material 
penetrates  the  temperature  jump  for  A0,-  =  1  or 
S^C.  However,  for  h  =  1000  m,  the  P  is  significantly 
reduced. 

A  more  realistic  temperature  distribution 
above  the  CBL  is  a  constant  dQa/dz.  Experiments 
simulating  this  distribution  as  well  as  a  jump  above 
a  well-mixed  layer  are  currently  underway  in  a  salt- 
stratified  tank  at  the  EP.4  Fluid  Modeling  Facility 
in  North  Carolina. 


W  (tons) 


Fig.  1.  Fraction  of  cloud  penetrating  an 
elevated  temperature  jump  (A0i,  °C)  as  a  function 
of  detonation  mass. 


Dispersion  parameters.  For  clouds,  Or  is 
dominated  by  entrainment  for  short  times  with 
(7r  =  arb  =  O.lSA/i.  At  intermediate  times 
{t  <  Ti),  the  CTr  may  be  dominated  by  ambient 
turbulence  in  the  inertial  subrange  with  err  ~ 
cTra  =  where  Ti  is  the  Lagrangian  time 

scale,  e  is  the  turbulent  kinetic  energy  dissipation 
rate,  and  ai  is  a  constant  (see  Thomson,  1990). 
At  long  times  {t  »  Tl),  <Tra  =  for 

homogeneous  isotropic  turbulence.  For  Ora,  we 
use  an  interpolation  expression  of  the  form  Cra  = 
f  {1  +  a2t /Tl)  to  satisfy  the  intermediate- 
and  long-time  results.  In  addition,  e  can  be 
written  as  e  =  ba^/Ti  in  homogeneous  isotropic 
turbulence. 

In  a  strong  CBL,  the  following  approximations 
can  be  made  for  z  >  O.lh:  e  ~  OAwl/h,  ~ 
0.6w,,  and  Tl  ~  0.7h/w,  (Weil,  1988).  These 
approximations  coupled  with  e  =  ba’^/Ti  lead  to 
b  =  0.78.  To  satisfy  the  long-time  dra  limit,  we 
must  have  02  =  0.62ai;  oi  is  estimated  to  be  0.57 
from  Thomson’s  two-particle  model  results.  The 
resulting  parameterization  for  Cra  in  the  CBL  is 


O'ra 

X 


0.36X3/2 

I'^i^o.six 


with 


_  w.x 

~  Hh 


(5) 


manner,  we  adopt  the  following  parameterization: 
tr3  =  For  clouds  dominated  by  buoyancy, 

^rb  — 

The  total  or  absolute  dispersion  is  necessary 
to  estimate  the  C  for  a  meandering  puff  or  plume. 
The  CTx  and  dy  in  Eq.  (2)  can  be  obtained  from  a 
parameterization  of  Taylor’s  theory:  d^  =  cr„t/(l-t- 
tl‘FrLxY/"^  and  similarly  for  dy.  The  Th  is  the 
Lagrangian  time  scale  for  the  u  component  and 
can  be  parameterized  by  T^x  oc  d^lh,  etc.  (e.g., 
see  Venkatram  and  Wyngaard,  1988).  For  the  CBL 
and  the  results  below,  we  use  Tix  =  '^Ly  —  0.7/i/rw, 
and  du  =  dr  =  O.Sra,. 


3.2  Some  Results 


We  have  computed  the  Cc  in  the  buoyant 
puff  and  the  mean  GLC  along  y  =  0  due  to  a 
meandering  puff  for  0.1  <  W  <  50  tons.  The  dr, 
dx,  and  dy  were  calculated  as  described  above.  In 
the  following,  the  cloud  buoyancy  is  characterized 
by  its  dimensionless  value 


Ft*  = 


Ft 

w'ihS  ' 


(6) 


we  used  w,  =2  m/s,  h  =  1000  m,  and  U  =  5  m/s. 

Figure  2  shows  the  dimensionless  concentration 
Cch^  /Q  as  a  function  of  X.  We  have  neglected 
cloud  penetration  of  the  inversion  but  included 
cloud  reflection  at  z  =  0,h  and  assumed  that 
he  =  Min  (Ah,  h).  The  large  vEiriation  in  the 
dimensionless  Cc  at  short  range  (X  <  1)  is 
due  to  the  buoyancy-induced  dispersion  drb-  As 
can  be  seen,  Cch^/Q  decreases  systematically  and 
significantly  with  an  increase  in  Ft,  due  to  the 
increase  in  drb  with  Ft,-  For  X  >  1,  the  curves 
converge  to  the  same  limit  because  at  long  times 
the  dr  is  dominated  by  dra,  which  is  independent 
of  Ft- 


X 


where  we  have  assumed  t  =  x/U. 

To  connect  the  short-,  intermediate-,  and  long-  Fig.  2.  Dimensionless  concentration  at  cloud 

time  relative  dispersion  regimes  in  a  continuous  centroid  versus  dimensionless  downwind  distance. 


Fig.  3.  SO2  concentration  at  cloud  centroid 
versus  downwind  distance. 


Fig.  4.  Dimensionless  mean  ground-level  con¬ 
centration  of  cloud  versus  dimensionless  downwind 
distance;  see  Fig.  3  for  key  to  lines. 


Figure  3  shows  dimensional  values  of  the  peak 
(Ce)  SO2  concentrations  in  the  cloud,  with  Q  = 
W  •  Ef  where  E/  {=  2.23  x  10“'*;  Andrulis,  1992) 
is  the  SO2  emission  factor.  In  Fig.  3,  the  order  of 
the  curves  is  reversed  from  Fig.  1 — the  curve  for 
W  =  50  tons  exhibits  the  highest  Cc-  The  reversal 
is  due  to  the  increase  in  Q  with  W,  which  overcomes 
the  decrease  in  Cc  due  to  the  increase  in  Crb  with 
Ft-  At  small  x,  all  of  the  curves  have  the  same 
slope:  Cc  cc.  because  Crb  oc  x*/^.  Some  curves 
exhibit  a  short  region  of  a  nearly  constant  Cc  with 
x;  this  is  due  to  puff  trapping  in  the  CBL.  At  large 
distances  (x  >  10  km),  clouds  for  all  cases  become 
uniformly  mixed  in  the  vertical  but  continue  to 
spread  lateradly;  thus,  Cc  oc  Q/cx^a  oc  Q/x  as 
shown. 

The  dimensionless  mean  GLC,  Ch^/Q,  along 
the  puff  centerline  is  shown  in  Fig.  4;  this 
mean  is  for  an  ensemble  of  meandering  puffs 
and  is  obtained  from  Eq.  (2)  with  reflection 
terms  included.  .A.gain,  the  highest  dimensionless 
concentration  occurs  for  the  smallest  Ft.',  this 
is  attributed  to  the  smaller  Ah  for  the  smaller 
detonations.  Likewise,  the  increase  in  the  distance 
to  the  maximum  concentration  with  Ft  is  due  to 
the  increase  in  Ah.  Note  that  for  X  <  1,  the 
Ch^JQ  can  be  two  orders  of  magnitude  smaller 
than  the  Cch^  fQ  (Fig.  2)  at  the  same  X  value,  but 
at  X  ~  10,  the  curves  from  both  figures  converge 
to  the  same  limit.  This  occurs  because  the  puff 
becomes  uniformly  mixed  in  z  and  the  CTi  ,  o-y  ~  Ora 
at  large  t  or  x. 

Figure  5  shows  the  mean  dimensional  GLC  for 
the  same  range  of  W  and  Ft  values  as  in  Figs.  2  - 
4.  Several  inter^^ng  features  are  found:  1)  A  non¬ 
monotonic  rarialion  occurs  in  the  maximum  GLC 
Cm  with  \V  and  Ft-  2)  The  variation  in  Cm  for 
0.1  <  W  <  50  tons  is  only  about  a  factor  of  4  even 
though  the  range  in  Q  is  a  factor  of  500;  the  weak 


dependence  on  Q  is  attributed  to  the  increase  in  Ah 
with  Ft-  3)  The  Cm  is  of  the  order  of  0.1  /rg/m^, 
which  is  the  lower  bound  for  Cc  in  Fig.  3. 

We  should  clarify  again  the  meaning  and  use 
of  C  in  Figs.  4  and  5.  It  is  the  meain  GLC 
along  y  =  0  due  to  an  ensemble  of  meandering 
puffs  and  probably  has  little  to  do  with  an 
observed  centerline  GLC  in  an  individual  puflf.  This 
computed  C  is  to  be  used  together  with  a  modeled 
cTc  in  a  concentration  probability  distribution  to 
estimate  the  peak  short-term  GLC  that  could  occur 
downstream  of  the  detonation.  The  peak  GLC 
would  correspond  to  some  specified  probability 
level. 

4.  SHORT-DURATION  RELEASES 


4.1  Dispersion  Model 

For  short-duration  releases  or  burns,  our  gen¬ 
eral  approach  is  an  integrated  puff  model  in  which 
the  short-term  mean  concentration  relative  to  the 
puff  centerline  is 


f*-  Qr  fit')  dt' 

0  (27r)3/2CT,.3,(T,.j,(7„ 


(7a) 


/  =  exp  - 


(76) 

where  t'  is  the  puff  emission  time,  tr  is  the 
total  release  duration,  Qr  is  the  continuous  source 
emission  rate,  z'  =  z  —  he,  drx  =  c^rx{t  —  t'), 
and  similarly  for  ary,cxrz-  The  integration  in  (7a) 
can  be  carried  out  analytically  for  limiting  forms 
of  —  t'),  etc.,  but  must  be  done  numerically 
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Fig.  5.  Mean  ground-level  SO2  concentration 
of  cloud  versus  downwind  distance;  see  Fig.  3  for 
key  to  lines. 

in  general.  Numerical  integration  is  required  when 
using  the  parameterization  =  <^ry  =  o'n  = 
a,  = 

The  integrated  puff  model  also  can  be  used 
for  the  mean  concentration  of  meandering  puffs  by 
replacing  the  relative  dispersion  by  the  absolute 
dispersion. 

In  the  following,  we  focus  on  the  Cc  for  a  short- 
duration  release  (burn)  and  consider  two  limiting 
cases.  1)  For  t  <  tr,  we  expect  the  rise  and 
dispersion  of  the  integrated  puff  to  reduce  to  that 
of  a  continuous  plume  for  sufficiently  strong  winds 
such  that  the  relative  dispersion  in  the  x  direction 
can  be  neglected.  2)  For  t  >  U,  the  C  field  should 
reduce  to  that  for  an  instantaneous  puff  but  with 
Q  =  Q^tr  and  Ft  =  {47T/Z)Fi,tr,  where  Ff,  (= 
Wor^gAQo/Q)  is  the  continuous  source  buoyancy 
flux.  As  will  be  shown  below,  the  Cc  for  the  long¬ 
time  puff  solution  is  lower  than  that  for  the  plume 
solution.  Thus,  we  take  the  plume  solution  as  an 
upper  bound  and  Cc  =  Min(Ccp/,  Ccjm),  where  Ccpi 
and  Ccpu  denote  the  Cc  values  for  the  plume  and 
puff,  respectively. 

The  mean  concentration  field  relative  to  the 
plume  centerline  is  given  by 

2-KUary‘ycz  ^  \  2cr2j,  2cr2^  ) 

Here,  the  plume  rise  is  attributed  to  buoyancy  and 
is  given  by  Ah  =  and  its  radius 

is  r  =  OAAh  (Briggs,  1984).  Source  momentum 
effects  can  be  included  in  the  future.  As  with  the 
puff  model,  we^^ll  assume  ayy  =  a-cz  =  o'r  and 
(Tr  =  “^ra  IS  given  by  Eq.  (5)  and 

the  plume  <Tr6  =  t/\/2  —  jU .  The 

Cc  —  from  Eq.  (8);  these  expressions 

can  be  expanded  to  include  reflection  at  2  =  0,  h. 


To  demonstrate  the  applicability  of  the  instan¬ 
taneous  puff  model  (Eq.  1)  for  long  times — t  >  tr 
cind  t  >  Ti,  we  carry  out  the  integration  in  Eq. 
(7a)  for  arx  =  (^ry  =  CTrz  =  CTr  =  (2<t2  gjid 

assume  =  (^w  We  ignore  the  dependence 

of  Ah  on  t' .  The  result  is 


Qr 


2alTL 


(x-r) 


where  erf  is  the  error  function  and  + 

+  {z  -  /ie)2.  We  evaluate  this  expression  at  a 
t  corresponding  to  the  center  of  the  cloud,  x  = 
U{t  —  tr/2),  OT  t  =  xjU  -t-  fr/2.  The  Cc  is  found 
to  be  Cc  =  (5rfr/[(27r)^'^^£Tr].  This  result  supports 
the  use  of  the  instcintaneous  puff  model,  with  Q  = 
Qrtr,  for  the  long-time  limit  of  a  finite-duration 
release. 

4.2  Some  Results 

Results  are  presented  for  the  dimensionless 
concentration  CcUh^fQr  for  the  plume  and  instan¬ 
taneous  puff  models,  with  reflection  at  z  =  0,  h 
included  in  both.  The  continuous  source  buoyancy 
flux  is  characterized  by  its  dimensionless  value: 


Figure  6  shows  the  dimensionless  Cc  for  the 
plume  model  with  F,  in  the  range  0.001  <  F.  < 
0.3.  The  trends  appear  similar  to  those  for  the 
puff  model  in  Fig.  2  although  the  variation  of 
CcUh? IQr  with  f,  is  not  as  great  as  for  the 
puff  model.  For  X  <  1,  the  decrease  in  the 
dimensionless  Cc  with  increasing  F.  is  due  to  the 
increase  in  Orb  with  Fb.  For  A'  >  1,  all  of  the 
curves  approach  the  same  asymptotic  curve;  this  is 
due  to  the  dominance  of  Ura  at  large  times  and  its 
independence  of  Fb. 

Figure  7  presents  the  dimensionless  Cc  for  both 
the  plume  and  puff  models  for  f,  =0.001  and  0.01 
amd  various  values  of  tr.  =  trW.jh.  The  time  scale 
hjw,  =  500  s  for  the  w,  (=  2  m/s)  and  h  (=  1000 
m)  used  here,  so  that  tr  ranges  from  50  s  to  500 
s  or  about  1  to  8  min.  The  plume  Cc  is  chosen 
as  long  as  it  exceeds  the  puff  Cc-  As  can  be  seen, 
the  distance  over  which  the  plume  solution  applies 
increases  as  tr.  does. 
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Fig.  6.  Dimensionless  concentration  at  plume 
centroid  versus  dimensionless  downwind  distance. 
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Fig.  7.  Dimensionless  concentration  at  plume 
and  puff  centroid  versus  dimensionless  downwind 
distance. 
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